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Abstract 

Propagation of light in a metamaterial medium which mimics curved spacetime 
and acts like a black hole is studied. We show that for a particular type of 
spacetimes and wave polarization, the time dilation appears as dielectric per¬ 
mittivity, while the spatial curvature manifests as magnetic permeability. The 
optical analogue to the relativistic Hamiltonian which determines the ray paths 
(null geodesics) in the anisotropic metamaterial is obtained. By applying the 
formalism to the Schwarzschild metric, we compare the ray paths with full-wave 
simulations in the equivalent optical medium. 

Keywords: Metamaterials, Anisotropic optical materials, Analogue gravity, 
Black hole, Schwarzschild spacetime 


1. Introduction 

One of the hot topics of modern technology is to build artificial materials 
whose permittivity and permeability can be properly engineered by incorporat¬ 
ing structural elements of subwavelength sizes. As a result, one can create ma¬ 
terials (called metamaterials) with the desired electromagnetic response which 
offers new opportunities for realizing such exotic phenomena as negative re¬ 
fraction, cloaking, super-lenses for subwavelength imaging, microantennas, etc. 
[!, 2 ]. 

Recently, it has also been recognized that metamaterials can be used to 
mimic general-relativity phenomena [3, 4]. The propagation of electromagnetic 
waves in curved spacetime is formally equivalent to the propagation in flat 
spacetime in a certain inhomogeneous anisotropic or bianisotropic medium [5- 
12]. Based on this equivalence, different general relativity phenomena have been 
discussed from the point of view of possible realization in metamaterials: optical 
analogues of black holes [13-18], Schwarzschild spacetime [12, 19], de Sitter 
spacetime [20-22], cosmic strings [23, 24], wormholes [25], Hawking radiation 
[26], the “Big Bang” and cosmological inflation [27, 28], colliding gravitational 
waves [29], among others. 

The deflection of light waves in gradient-index optical materials mimicking 
optical black holes was studied theoretically [13, 14, 18] and experimentally 
[15-18]. These materials, called by authors “omnidirectional electromagnetic 
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absorbers”, are characterized by an isotropic effective refractive index. A real 
cosmological black hole (BH) can often be described by an anisotropic spacetime, 
as, for example, the case of the Schwarzschild BH. For that case, one should 
determine the permittivity and permeability tensors instead of the refractive 
index in order to introduce the equivalent optical medium [12, 19]. Chen et al. 
[19] simulated the wave propagation outside the Schwarzschild BH and observed 
in their numerical results the phenomenon of “photon sphere”, which is an 
important feature of the BH system. It would be interesting to go further 
and study the propagation of light waves in optically anisotropic media which 
mimic cosmological BHs and compare the results with ray paths obtained from 
the Hamiltonian method. 

The aim of this letter is twofold. First, we determine the constitutive rela¬ 
tions of an inhomogeneous anisotropic medium which is formally equivalent to 
the static spacetime metric obeying rotational symmetries and can be applied, 
in principle, to the medium either in isotropic or anisotropic form. Second, by 
making use of the eikonal approximation to the wave equation, we obtain the 
expression for the optical Hamiltonian which we found to be identical to the one 
obtained from general relativity for null geodesics, but different from the opti¬ 
cal Hamiltonian used in Refs. [13, 15, 30, 31]. Then we apply the formalism to 
the Schwarzschild spacetime that is a solution to the Einstein held equations in 
vacuum [32] . We compare the wave propagation with the ray dynamics outside 
the BH in the effective medium and obtain a very good correspondence. As an 
interesting feature we find that light does not propagate in the direction of the 
wave normal, there is an angle between the wave velocity and the ray veloc¬ 
ity. The obtained results are discussed from the point of view of metamaterial 
implementation. 

2. General relativity in a metamaterial medium 

2.1. Medium parameters 

Long time ago, Tamm pointed out the parallels between anisotropic crystals 
and curved spacetimes [5]. Later studies showed [6-11] that the propagation of 
light in empty curved space distorted by a gravitational held is formally equiva¬ 
lent to light propagation in hat space filled with an inhomogeneous anisotropic 
medium. 

Indeed, consider a spacetime background with a general metric 1 

ds 2 = goo At 2 + 2g 0i d t Ax 1 + gtj Ax 1 dad, (1) 

where i,j = 1,2,3 run over arbitrary spatial coordinates. Then, it can be shown 
[11] that the covariant Maxwell’s equations written in curved coordinates can 
be transformed into their standard form for hat space but in the presence of 


1 From now on we follow the standard notations for covariant (subindices) and contravariant 

(superindices) quantities. 
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an effective medium. The constitutive relations of the equivalent medium have 
been found in the form [7]: 

D i = £ ij Ej - (T x H)*, B i = fjiHj + (fx E)*, (2) 


which connect the fields D, B, E and H via nontrivial permittivity and perme¬ 
ability tensors 


£ lJ = n' 3 = 


5oo 


( 3 ) 


and a vector T given by 



( 4 ) 


Here, is the inverse of and g is the determinant of the full spacetime metric 
with /x, v = 0,1, 2, 3. Note that the information about the gravitational field 
is essentially embedded in the material properties of the effective medium: the 
tensors £ lJ . p lJ which are symmetric and should be equal, and the vector T 
which couples the electric and magnetic fields. The invention of metamaterials 
during the last decade [1, 2] opened up the possibility to design electromagnetic 
media corresponding to different spacetimes [3, 4, 19-29]. 

In this letter, we consider a static spacetime metric associated with a spher¬ 
ically symmetric cosmological BH. Due to time-reversal symmetry, got = 0 and 
the coupling between the electric and magnetic fields vanishes, T = 0. The 
metric (1) in (t, r, 9, ip) coordinates can then be written in a generic form as [32] 


ds 2 = goo(r) d t 2 + g rr {r) {dr 2 + /(r) [r 2 d 9 2 + r 2 sin 2 9dp 2 ]] , (5) 

where / is the “anisotropic factor”. Note that the metric (5) obeys the rotational 
symmetries in the three-dimensional (r,6,ip) space. 

Then, we have to project the metric (5) into a flat background to obtain 
the medium parameters in the Cartesian coordinate system. To do that, we 
apply a coordinate transformation and, from Eq. (3), we get the permittivity 
and permeability tensors in the form: 


/ 9rr 

rpi rpj 

V 900 



( 6 ) 


where is the Kronecker delta, r = x 2 + y 2 + z 2 , and we denoted the Carte¬ 
sian coordinates (a: 1 ,^ 2 ,^ 3 ) = ( x,y,z ). It is seen that whenever f ^ 1, the 
permittivity and permeability tensors contain the off-diagonal elements and the 
equivalent medium is essentially anisotropic. Only in the case of / = 1, the off- 
diagonal elements vanish and the medium becomes completely isotropic with all 
the diagonal elements equal to the refractive index: n(r) = \J —g r r/ goo- Note 
that in general relativity the spacetime with / = 1 in Eq. (5) is said to be con¬ 
formal to flat space. Every static spherically symmetric spacetime with / ^ 1 
can be converted to conformally flat form by an appropriate transformation of 
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the radial coordinate: r —> p. The new radial coordinate is obtained by 


p = r exp 



Vf( r ') 



( 7 ) 


where the isotropic boundary condition at infinity, /(oo) = 1, is taken into 
account. The line element in the (t,p,6,<p) isotropic coordinates takes the con¬ 
formally flat form: 

ds 2 = gooK/o)] df 2 + A(p) ( dp 2 + p 2 sin 2 9 dip 2 + p 2 d9 2 ), (8) 

where the time dilation term g 00 and the conformal factor A = g rr fr 2 /p 2 are 
calculated by means of the function r(p) which should be obtained by inverting 
(7). Thus, the permittivity and permeability tensors are simply reduced to the 
isotropic refractive index: 


e ij = p ij =S ij J-— = S ij n(p ). (9) 

V 9oo 

The equivalent medium determined by Eq. (9) is still inhomogeneous since the 
refractive index varies with the radial coordinate, but the light velocity in the 
medium becomes isotropic, a property that is much simpler to implement in 
metamaterial design. 

In what follows, we will compare the results for light propagation in two 
different equivalent media - isotropic and anisotropic - both corresponding to 
the same spacetime metric in order to see the physical differences. 

2.2. Electromagnetic fields. TE and TM polarizations 

The results from the previous section indicate that an electromagnetic field 
can be thought of as propagating in flat background but in the presence of a 
medium whose properties are constructed from a curved spacetime. The fields, 
for the static case we consider, are related by: 

D l =£ ij Ej, B' = pP H 3 . (10) 

Due to the anisotropy of the medium for the metric in nonconformally flat form, 
the electric displacement field D is not in the direction of E, and the magnetic 
induction field B is not in the direction of H. To simplify the treatment of the 
problem, we consider the propagation of light in the equatorial plane, z = 0. In 
such a case, one of the anisotropies - electric or magnetic - can be eliminated. 

Indeed, consider the TE polarization for an electromagnetic wave for which E 
is perpendicular to the x-y plane. Equation (6) for z = 0 leads to £ xz = £ yz = 0, 
hence the directions of D and E coincide. This means that £ zz is the only 
relevant matrix element which connects the nonzero electric components of the 
field: D z = £ ZZ E Z , and the electric anisotropy of the medium is irrelevant. As 
for the magnetic components, we obtain p xz = p yz = 0, B z = H z = 0, and 
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B x = p xx H x + l-i xy H y and B y = p xy H x + n vv H y . Hence, B and H are confined 
to the x-y plane, but their directions do not coincide. 

Similarly, one can consider the TM polarization for which H is perpendicular 
to the x-y plane. In such a case, one gets B z = p zz H zl D x = e xx E x +£ xy E y , and 
D v — s xy E x +£ vv E y . This means that the directions of B and H coincide, and 
the magnetic anisotropy is irrelevant. It should be noted that, since £ zy = 
and they are both symmetrical, the TE and TM polarizations are completely 
equivalent. These two cases are summarized in Fig.l. 



(a) Magnetic anisotropy. (b) Electric anisotropy. 


Figure 1: Electromagnetic fields H, B, E, D, wavevector k and Poynting vector S for: 
(a) TE polarization in a magnetically anisotropic medium, and (b) TM polarization 
in an electrically anisotropic medium. 

The medium anisotropy has very important physical consequences on the 
propagation of the electromagnetic waves. If we look at the monochromatic 
plane wave with wavevector k, then D, B and k form an orthogonal vector 
triplet. On the other hand, E, H and the Poynting vector S form another 
orthogonal vector triplet. Therefore, the energy does not propagate in the di¬ 
rection of the wave normal [33] and the ray velocity and the wave velocity are 
not equal, unlike what happens for an isotropic medium. It is seen that the 
angle between S and k (denoted by a in Fig.l) is equal to that between B and 
H in the TE polarization or between D and E in the TM polarization. For the 
angle a we find: 


p xy __£ xy _ (1 - f)xy 
p yy £yy x 2 + fy 2 


( 11 ) 


Only when / = 1, the medium becomes isotropic, a vanishes and the ray and 
wave velocities become equal. 


2.3. Schwarzschild spacetime 

Equations (6) and (9) can now be applied to any static spherically symmetric 
spacetime. As an example, we consider the Schwarzschild metric that is a 
solution to the Einstein field equations in vacuum [32] . This metric in spherical 
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coordinates reads: 


ds 2 = — ^1--^ dt 2 + ^ + r 2 (sin 2 9 dip 2 + dd 2 ). (12) 

r 

Here, the speed of light c = 1, the Schwarzschilcl radius r s = 2 GM/c 2 , G is the 
universal gravitational constant and M is the mass of the source of gravitation. 
Comparing with Eq. (5), one can easily identify: g 00 = —(1 — r s /r), g rr = 
— l/ffooj and 1 — / = r s /r. Since the Schwarzschild metric (12) is not in a 
conformally flat form, one obtains the anisotropic equivalent medium tensors 


e tJ = ^ 




(13) 


and from Eq. (11) we find the angle between the wave and ray velocities: 


tana = 


xy r s 

r 3 _ y2 Is 


(14) 


Equation (13) is a compact version of the medium tensors which can also be 
found as matrices in Refs. [19, 22]. 

The Schwarzschild metric can be written in a conformally flat form by means 
of the well known substitution (see, e.g. [32]) 


r = p 



The inverted formula for the isotropic radial coordinate p reads: 



(15) 


(16) 


which can also be obtained from Eq. (7) by substituting / = 1 — r s /r. In 
isotropic coordinates, the Schwarzschild metric becomes [32] 


ds 2 


dt 2 +(l + ^j\dp 2 + p 2 sm 2 6dp 2 + p 2 de 2 ). (17) 


Finally, the refractive index of the equivalent medium is found from Eq. (9) as 
(see also Ref. [9]) 



(18) 


The formulas given in this section for the Schwarzschild spacetime will be used 
in Sec. 5 for numerical simulation. 
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3. Maxwell’s equations in anisotropic medium 


The starting point of our analysis of wave propagation in the equivalent 
medium is the Maxwell equations written in a source-free form 

r) 1~) 

V • B = 0, V-D = 0, V x E = ——, V x H = . (19) 

dt' dt v 

These equations should be supplemented by the constitutive relations (10) with 
the medium tensors (6) expressed through the spacetime metric coefficients. If 
the operating frequency bandwidth is sufficiently narrow, the dispersion can be 
neglected and one can consider a monochromatic wave with frequency ui. Equa¬ 
tions (19) can then be reduced [34] to the wave equation for a time-harmonic 
electric field vector 


V x 



(V x E u ) 




( 20 ) 


The medium is generally anisotropic, hence, e and /.i 1 denote the permittivity 
and inverse permeability tensors, respectively. 

As it was discussed in Sec. 2.2, the electromagnetic propagation in the equa¬ 
torial plane, z = 0, can be described in terms of either TE or TM waves. For 
the TE wave, E u (x,y,z) = E(x,y)e~' kzZ z, we obtain the equation for the 2 - 
component of the electric field E in the two-dimensional ( x , y) space 


d 

f dE 

dE\ 

d 

( dE 

dE\ 

dx 

\ xv ~dy 

~ Myy d^J 

dy 1 

f^XX~P{ 

V d v 

~^ y ~di) 


U1 2 £ ZZ E, 


( 21 ) 


where fi xx , Hxy, Hyy are the corresponding components of the tensor p 1 and 
e zz corresponds to e. Since only one component of e has entered into the final 
equation (21), the permittivity can be chosen isotropic with all its diagonal 
elements equal to £ zz and the off-diagonal terms equal to zero. In this case, the 
equivalent medium will be magnetically anisotropic [see Fig. 1(a)]. 

The wave equation for the TM case can also be easily obtained from the 
above formulas by means of the substitutions: /W, E —► H. In such a 

case, the medium is electrically anisotropic [Fig. 1(b)], and one should solve the 
equation for the magnetic field H(x, y). 

We would like to point out that we do not use the plane wave approximation 
to treat the problem, instead, we will solve the wave equation (21) numerically 
for a Gaussian beam with the appropriate boundary conditions. This situation, 
as we believe, is more realistic from the experimental point of view. 


4. Hamiltonian formulation 

In the limit of geometrical optics, for which the wavelength is much smaller 
than the scale of variation of the medium parameters, electromagnetic waves 
should follow ray trajectories (null geodesics of the background spacetime). To 
test the validity of our scheme, it would be interesting to compare the full-wave 
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numerical calculation of the propagation of light with its geometrical optics 
limit, for which analytical solutions exist. To accomplish this task, we have to 
define the Hamiltonian. 

Let us start with the wave equation (20) , which we apply again to the equa¬ 
torial plane, z = 0, and consider the TE polarization for the wave propagation. 
Given that the medium tensors are diagonal in the spherical coordinate system, 
it is advantageous to rewrite Eq. (21) in polar coordinates: E(x,y ) —> E{r,<p). 
One gets, 

J_d_ 
r 2 dp 

where we denote p r = p rr , = // w , and e z = £ zz , which are the medium 

parameters in polar coordinates. 

In the eikonal approximation we seek the solution in the form: E(r,p) = 
E 0 exp[«5'(r, </>)], where E 0 is the slowly varying amplitude and S(r, ip) is the fast 
varying eikonal [35]. Substituting into Eq. (22) and neglecting the small terms 
corresponding to the derivatives over slowly varying functions, we obtain: 


1 9E\ 2 „ „ 

- T, — ) + Ul £ z E — 0. 

Pr Oip 


( 22 ) 


1 d ( r dE 
r dr \p, v dr 


1 

Pip 


(dS\ 

I I 1 1 

(d S \ 

\dr ) 

p r r 2 

{dtp] 


£ z OJ 2 = 0 , 


(23) 


which can be rewritten in terms of the spacetime metric coefficients as 


1 

Qrr 



fr 2 \dp) 


— cc 2 = 0. 

goo 


(24) 


But this is nothing else than the eikonal equation for a light ray in a gravitational 
field [35] 
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ot(3 


dS dS 


= 0 


(25) 


dx a dx& 

with the indices a, /? running over space and time. The Hamiltonian representa¬ 
tion of Eq. (25) is obtained by rewriting it in terms of the variables p a = dS/dx a , 
where p a is the 4-momentum. One gets 


H = g a/3 p a pp = P^Pp = o. 


(26) 


Thus, with the definitions p r = dS/dr , p v = dS/dip , and pt = u, just plugging 
them into Eq. (24), we obtain 


H = 


1 9 1 

Pt ^ 

$00 9rr 




(27) 


The optic Hamiltonian H given by Eq. (27) describes the ray trajectories in the 
equivalent anisotropic medium and is equivalent to the one used to determine the 
null geodesics in the gravitational field. It takes into account the time dilation 
by means of the term goo, the spatial curvature by the metric coefficient g rr and 
the medium anisotropy through the function /. The latter should result in the 
deviation of the momentum vector p from the ray propagation direction. 






Having found the Hamiltonian T~L , the ray dynamics can easily be calculated 
from Hamilton’s equations of motion: 


d'H . dH 

~ dp a ’ Pa ~ dq a 


(28) 


where the canonical coordinates in our case are: q a = p a = ( pt,Pr,Pip ) 

and the “dot” denotes the derivative over the affine parameter which varies 
along the trajectory. Note that in the Hamiltonian (27), the coefficients goo, g rr 
and / are functions of only the radial coordinate. Therefore, there exist two 
integrals of motion, p t = 0 and p v = 0, which just state the conservation of 
energy and angular momentum for the static spherically symmetric system. In 
principle, these two constants of motion can be introduced; but in fact, since 
photons are massless, the ray trajectories are determined by only one constant 
- the impact parameter b which is the ratio between them [36]. 

We can now determine the ray trajectories for a particular spacetime metric 
by means of Eqs. (28) with the appropriate boundary conditions. To be specific, 
we consider the Schwarzschild metric (12) for which the Hamiltonian can be 
readily obtained by means of Eq. (27) 


^ — ~ l — is. (l 

r 




2 


(29) 


On the other hand, in isotropic coordinates (see Eq. (17)), the Schwarzschild 
Hamiltonian becomes 


H = - 





(30) 


For both cases, isotropic and anisotropic, we will compare different kinds of tra¬ 
jectories with the full-wave numerical simulations in Sec. 5. Note that Eq. (26), 
when applied to the optical medium, gives the generalized dispersion relation 
[4] also known as the Fresnel equation [33, 34]. The discussion of this relation 
for the Schwarzschild equivalent medium can be found, e.g., in Ref. [12]. 

It would be appropriate to make some important remarks here. 

(i) By comparing the terms in Eqs. (23) and (24), one can see that for the TE 
wave, the time dilation term g 0 o is determined by the dielectric permit¬ 
tivity of the equivalent medium, while the spatial metric components g rr 
and / depend on the magnetic permeability. In other words, the time and 
space metric coefficients correspond to different properties of the medium: 
either electric or magnetic. Similarly, if we consider the TM wave, one 
should replace: e l ' J // J ’; E —>• H (see Sec. 2.2), and the correspondence 
will be just the opposite: goo is related to p z , while g rr and / correspond 
to e r and e v . 

(ii) Equation (23) can be divided by e z and the medium parameters can be re¬ 
defined according to: jl v = p^,£ z , p, r = p r £z and e = 1. It is easy to check 
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that the new medium gives exactly the same ray trajectories, although it 
has the permittivity e of free space and the inhomogeneous permeability 
fi. Actually, the Hamiltonian T~L may be multiplied by any arbitrary func¬ 
tion of the coordinates (by goo? i n this case) without changing the ray path 
obtained from the equations of motion since only the parametrization of 
the affine parameter changes. This means that in the geometrical optics 
limit, one can introduce a set of equivalent media, all giving the same light 
ray paths. But the wave equation (22) (or in general Eq. (21)) cannot be 
renormalized in this way, and therefore the interference patterns for those 
analogous media will be different. 


5. Comparison of wave propagation with ray dynamics 

In this section we present the results obtained numerically for the propaga¬ 
tion of light waves in an anisotropic medium which mimics a static spherically 
symmetric cosmological spacetime. As an example, we consider the Schwarzschild 
spacetime with the goal to compare the results for two sets of coordinates: spher¬ 
ical and isotropic. On the other hand, to validate our theoretical approach, we 
wish to compare the full-wave numerical simulation with the ray trajectories 
obtained within the Hamiltonian framework. 

We solve Eq. (21) in a rectangular 2D geometry of (x,y) space for a TE 
polarized wave injected from the right (see Figs. 2-3). The medium parameters 
we use are given by Eq. (13) for the anisotropic spacetime and by Eq. (18) for 
the isotropic spacetime. In both cases, we solve the equations in Cartesian 
coordinates. The computational domain is surrounded by a perfectly matched 
layer that absorbs the outward waves to ensure that there are no unwanted 
reflections, and the simulations are done by means of a standard software solver. 

It is known that the Schwarzschild black hole presents an event horizon at 
r = r s (in isotropic coordinates, at p = r s /4 = p s ) where the gravity is so 
strong, that light cannot escape [36]. The metric has a coordinate singularity 
there, which translates into singular medium parameters. To avoid it, we set an 
effective horizon at r s + S with 5 being small positive number 2 and impose an 
absorbing inner core for r < r s + 5 [13, 19]. 

We use a Gaussian shape for the TE wave injected at the boundary, however, 
care should be taken for the injection direction. If the medium is anisotropic, the 
beam does not propagate in the direction of the wave normal, as we discussed 
in Sec. 2.2. There will be an angle a between the wave velocity and the ray 
velocity (see Fig. 1). To ensure that the wave is always injected in a controlled 
way, in our case along the cc-axis, independently of the point of injection (that 
is, for different impact parameters 6), we impose: n = — (cosa)x — (sina)y for 
the normal unit vector n of the wave front. The angle a is determined by Eq. 


2 We set 5 = 0.05 r s (anisotropic case) and 8 = 0.11 r s (isotropic) in the calculations. We 
have checked that the variation of <5, whenever it is small enough, does not affect the global 
distribution of fields in the whole domain except some tiny shell at the horizon. 
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(11) applied to the point of injection (xo,yo). The fields correspond to those 
depicted in Fig. 1(a), in particular, the electric field vector is perpendicular to 
the plane of simulation and the medium has a magnetic anisotropy. Figure 1(b) 
corresponds to the case of TM wave injection into an electrically anisotropic 
medium. This case, in principle, can also be simulated. Because the injected 
wave is not a plane wave, the wavevector k in Fig. 1 should be replaced by the 
wave normal vector n of the Gaussian beam. 

The impact parameter b is the key quantity which allows us to distinguish 
between different types of ray trajectories [36]. Depending on its value, one may 
observe capture (6 < b c ) or deflection (b > b c ). Its critical value b c = 3-^/3r s /2 
determines the so-called photonic sphere, which corresponds to an unstable 
circular orbit of radius 3 r s /2. To relate the value of b with the initial coordinates 
(x 0 ,yo), we use the following formula 


b = 2/o 


rof 


9oo y/xl + fy\ 


(31) 


with the metric parameters taken at tq = \Jx\ + j/q . Equation (31) has been 
derived under the assumption y = 0 at the injection point. It can be verified 
that at infinity (tq —1 oo) b approaches yo as it should according to the definition 
of the impact parameter. In anisotropic spacetimes, the injection momentum 
p = p x x + p v y is not collinear with the propagation direction. Since we impose 
that initially the ray propagates along the £-axis, at the injection point the 
momentum forms an angle given by 


PyQ 
Px 0 


(1 ~ f)x 0 y 0 
x o + fy o 


tan a. 


(32) 


It is exactly the same angle a between the wave normal and the Poynting vector 
and between the B and H fields, all given by (11). Here we clearly see how the 
spacetime anisotropy results in the anisotropy of the medium. 

In Fig. 2 we show the results of numerical simulations for Gaussian beams 
incident with different impact parameters b for the case of a Schwarzschild 
anisotropic spacetime. The ray paths calculated from the Hamilton equations 
are superimposed in white on the wave pattern (shown in color) simulated by 
solving the wave equation. We observe a good correspondence between the two 
approaches: the ray path follows the center of the Gaussian beam even when 
it deflects. Moreover, the ray trajectories are consistent with the behavior ex¬ 
pected from theory: we observe the ray capture for b < b c [Figs. 2(a), (b)], and 
its deflection for b > b c [Figs. 2(c), (d)]. The wave solutions present a more 
complex behavior due to interference effects. The Gaussian beam splits into a 
set of rays or “subbeams”, one part of it is captured due to having b < b c and 
another part bends around the photonic sphere of the black hole and interferes 
with the primary beam. This effect is more pronounced in Figs. 2(b), (c) for 
the values of b which are not very far from the critical value b c « 2.598r s . It is 
interesting to see that the subbeams are regularly spaced while bending around 
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the black hole and some of them exhibit smaller subwavelength ripples. Similar 
behavior was observed in Ref. [19] for a larger value of the wavelength. 

In Fig. 3 we present the results for the isotropic case of the Schwarzschild 
metric using the same set of the impact parameters b as for the anisotropic 
case of Fig. 2. Again, we observe a good agreement between the wave prop¬ 
agation and the ray tracing. The principal difference is that now the local 
velocity of light becomes isotropic and it is determined by a scalar function - 
the refractive index n(p). As a consequence, the interference pattern we ob¬ 
serve is substantially different. One can see irregular “jets” leaving the system 
in radial directions, similar to those obtained by Genov et al. for the isotropic 
spacetimes [14]. The small ripples are also clearly seen. Since we simulate an 
isotropic medium, another distinctive feature we observe, is that the wavefronts 
are perpendicular to the beam propagation, i.e. a = 0. 

On the other hand, if we compare the ray trajectories for the isotropic and 
the anisotropic cases, they are rather similar: the same kinds of orbits (capture 
or deflection) are observed corresponding to the same value of the impact pa¬ 
rameter b. Notice that the parameter b is a coordinate-independent quantity, 
since it is defined as a ratio between the angular momentum and energy of the 
photon measured at infinity. Other physical properties, like the deflection angle 
between the incoming and outgoing light rays measured by a distant observer 
at infinity should also be equal independently of the coordinate system used. 
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Figure 2: TE Gaussian beam of wavelength A = 0.15r a compared with ray path (super¬ 
imposed in white) in metamaterial mimicking the Schwarzschild anisotropic spacetime 
for different impact parameters: (a) b = r s , (b) b = 2r a , (c) b ~ b c (slightly above), and 
(d) b — 3r a . The event horizon r a and the photon sphere of radius 3r a /2 are depicted 
by black circles. 
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Figure 3: TE Gaussian beam of wavelength A = 0.15r a compared with ray path (super¬ 
imposed in white) in metamaterial mimicking the Schwarzschild isotropic spacetime 
for different impact parameters: (a) b = r s , (b) b = 2r s , (c) b ~ b c (slightly above), 
and (d) b = 3r a . The event horizon p a and the photon sphere of radius (2 + y / 3)r a /4 
are depicted by black circles. 
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6. Concluding remarks 

In conclusion, we have derived the constitutive relations of an inhomogeneous 
anisotropic medium which is formally equivalent to a class of static spacetime 
metrics associated with a spherically symmetric cosmological black hole. It is 
known that every spacetime metric can be written in different coordinate sys¬ 
tems and therefore, can be projected in many different ways into the correspond¬ 
ing effective medium, each with specific material properties. This procedure is 
at the basis of the transformation optics [3, 4]. In particular, a static spheri¬ 
cally symmetric spacetime written in arbitrary (non-isotropic) coordinates can 
be equivalently described in isotropic coordinates (conformally flat form) after 
performing the appropriate coordinate transformation. 

We have analyzed the medium properties and the propagation of light through 
the effective media corresponding to these two cases of interest: anisotropic and 
isotropic. For the anisotropic spacetime, we found that only one kind of the 
medium anisotropy would be essential: either magnetic or electric, for TE or 
TM polarized waves respectively. In both cases, light does not propagate in the 
direction of the wave normal. There appears an angle between the wave velocity 
and the ray velocity, specified by Eq. (11), which is related to the anisotropic 
factor of the metric. For some spacetimes, this angle can even be negative giv¬ 
ing such interesting phenomenon as negatively refracting medium [37]. However, 
this cannot be observed for the metric (5). 

It is interesting to note that for the kind of metric we consider, the time and 
space metric coefficients correspond to different properties of the medium: either 
electric or magnetic. For instance, for the TE wave, the time dilation term de¬ 
termines the dielectric permittivity of the equivalent medium, while the spatial 
metric components correspond to its magnetic permeability. This phenomenon 
cannot be seen in the isotropic medium characterized by the scalar refractive 
index. We presented the results for the TE wave in the magnetically anisotropic 
media, although the case of the TM wave in the electrically anisotropic media 
perhaps would be easier to realize in the laboratory experiment. Both consid¬ 
erations are theoretically equivalent. 

By applying the eikonal approximation to the wave equation, we have also 
obtained the expression for the optical Hamiltonian which was used to simulate 
the ray paths. It coincides with the Hamiltonian obtained from general rela¬ 
tivity for null-geodesics. It is clearly seen from our theoretical approach, that 
despite the fact that in the geometrical optics limit different media may give 
the same ray paths, their full-wave description will differ. Taking as an example 
the Schwarzschild spacetime metric, we obtained a very good correspondence 
between the wave propagation and the ray trajectories for the anisotropic and 
isotropic media. Although the isotropic case is easier to implement in meta¬ 
materials, we believe that due to a rapid progress in metamaterial technology 
[38] (see also proposal [22]), anisotropic media will also be available to model 
the cosmological phenomena in the laboratory and many physically interesting 
phenomena could be observed. 
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